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Abstract. We present a determination of the optical/UV AGN luminosity function and its evolution, based on a large sample 
of faint (R < 24) QSOs identified in the COMBO-17 survey. Using multi-band photometry in 17 filters within 350nm & ?i (, s 
930 nm, we could simultaneously determine photometric redshifts with an accuracy of a z < 0.03 and obtain spectral energy 
distributions. The redshift range covered by the sample is 1.2 < z < 4.8, which implies that even at z — 3, the sample reaches 
below luminosities corresponding to Mb = —23, conventionally employed to distinguish between Seyfert galaxies and quasars. 
We clearly detect a broad plateau-like maximum of quasar activity around z — 2 and map out the smooth turnover between 
z — 1 and z — 4. The shape of the LF is characterised by some mild curvature, but no sharp 'break' is present within the range 
of luminosities covered. Using only the COMBO-17 data, the evolving LF can be adequately described by either a pure density 
evolution (PDE) or a pure luminosity evolution (PLE) model. However, the absence of a strong L*-like feature in the shape 
of the LF inhibits a robust distinction between these modes. We present a robust estimate for the integrated UV luminosity 
generation by AGN as a function of redshift. We find that the LF continues to rise even at the lowest luminosities probed by our 
survey, but that the slope is sufficiently shallow that the contribution of low-luminosity AGN to the UV luminosity density is 
negligible. Although our sample reaches much fainter flux levels than previous data sets, our results on space densities and LF 
slopes are completely consistent with extrapolations from recent major surveys such as SDSS and 2QZ. 

Ksy words, surveys — galaxies: active — galaxies: Seyfert — quasars: general 



1. Introduction 

The luminosity function of quasi-stellar objects (QSOs) and 
its evolution with redshift provides one of the most im- 
portant tools for the cosmic demography of active galac- 
tic nuclei (AGN). It constrains physical models for QSOs, 
particularly those for the growth of supermassive black 
holes in galaxies within the conte xt of hierachical col- 
lapse of structure in the un iverse (lHaehnelt & Rees 19931 
lHaiman & Loeb 1998llKauffmann & Haehnelt 2000l) . It is also 
relevant for un derstandi ng the extragalactic UV background 
jMeiksin & Madau 1993llBovle & Terlevich 19981) . 

Previous studies of the QSO luminosity function (QLF) 
have established a strong rise in the activity with look back tim e 
from the local universe to redshifts z ~ 2 jBovle etal. 1 988). 
The main debate at these low to intermediate redshifts is now 
about the question whether the shape of the QLF changes with 
redshift JBovle et al. 200ri Iwisotzki 20001 iMivaii et al. 2 000: 
ICowie et al. 20031) . If so. this would mean that low-luminosity 
AGN evolve differently from high-luminosity QSOs, possibly 



calling for a substantial revision of our understanding of the 
cosmic duty cycle of nuclear activity in galaxies. 

At higher redshift z > 2, the obvious expectation that the 
rise of the QSO activity with redshift mus t turn over at some 
point was satisfied b y observations at z > 3 JWarren et al. 19941 
ISchmidt et al. 19951 hereafter WHO and SSG). However, opti- 
cal surveys in this redshift regime were often plagued by se- 
lection effects an d small number sta tistics. Moreover, studies 
of X-ray-selected jMivaii et al. 20 00) and radio-selected QSOs 
( Jar vis & Rawlings 200fj|) did not strongly support claims of a 
very steep drop in nuclear activity for z > 3, keeping the issue 
unresolved whether the detected turnover was physical reality 
or not. A rather robust observation of the space density decline 
beyond z > 3.6 has been established recently, albeit only for 
very luminous QSOs, by the Sloan Digital Sky Survey (SDSS, 
iFan et al. 200lh . 

Most optically selected QSO samples are expected to be 
more or less complete at either z 2.2 or z 3.6, where 
QSOs show conspicuous colours in broad-band searches, while 
within this redshift band, confusion arises with stars and com- 
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pact low-redshift galaxies. Any practical approach of following 
up QSO samples in this intermediate redshift range is forced to 
avoid strong contamination and observed mostly extreme ob- 
jects, resulting in a high degree of incompleteness that is diffi- 
cult to account for. Thus, while the fact that cosmic QSO activ- 
ity shows a maximum around z — 2-3 is not doubted as such, 
this maximum has very rarely been detected in a single survey. 

This selection issue was one of the original motivations for 
the COMBO-17 survey. This survey is based on photometric 
classification and redshift estimation using a set of 17 filters, 
of which 12 are medium-band filters with ~ 10,000 km/sec 
FWHM, well matched to optimum detection of QSO emission 
lines and spread across the range of wavelengths observable by 
modern CCD detectors. 

Another major driver for this project was to probe fainter 
regimes of the AGN luminosity function, a task which has al- 
ways beed limited by the overwhelming need for spectroscopic 
follow-up telescope time. Using medium-band spectrophotom- 
etry, a reasonably complete and clean AGN sample could be 
obtained across a wide range of redshifts down to R ^ 24. 
Based on such a deep sample, we were able to derive luminos- 
ity functions down to Mb —23 at redshifts above 3, entering 
the domain of Seyfert galaxies even at high redshift. 

In this paper, we present a new determination of the lu- 
minosity function of faint optically selected QSOs, aimed at 
a broad redshift range around the elusive turnover. The anal- 
ysis is based on a sample of 192 objects between z — 1.2 and 
z — 4.8, all selected from the COMBO-17 survey. The survey is 
briefly described in Sect. 2, followed by a discussion of techni- 
cal aspects such as completeness and redshift quality in Sect. 3. 
In Sections 4 through 7 we present and discuss the results. 

2. COMBO-17 observations 

The COMBO-17 project ("Classifying Objects by Medium- 
Band Observations in 17 Filters") was designed to provide a 
sample of ^50,000 galaxies and several hundred AGN with 
precise photometric redshifts and spectral energy distributions 
(SEDs). As shown below, the filter set provides a redshift ac- 
curacy of g 7 ~ .03 for quasars (and similarly for galaxies; cf. 
Iwblf et al. 20031) . smoothing the unknown true redshift distri- 
bution of the sample only very mildly and certainly allowing 
the derivation of luminosity functions. 

The first step of the COMBO-17 data analysis was to con- 
vert the photometric observations into a very low resolution 
'fuzzy spectrum', allowing for simultaneous spectral classifi- 
cation into stars, galaxies and QSOs, as well as for accurate 
redshift and SED estimation for the latter two categories. The 
full survey catalogue will contain about 50,000 objects with 
classifications and redshifts covering a solid angle of 1.0 deg 2 . 
This fuzzy spectroscopy consciously compromises on redshift 
accuracy in order to obtain large samples of quasars with a rea- 
sonable observational effort. 

While the photometric redshift technique has already been 
appl ied to gal axy samples about 40 years ago jBaum 1 962: 
Butc hins 1983b . we have modified and improved the approach 
by increasing the number of filters and narrowing their band- 
width to obtain better spectral resolution and more spec- 



Table 1. The COMBO-17 filter set: Exposure times and 10c 
magnitude limits reached for point sources, averaged over all 
three fields. The 7?-band observations were selected to be taken 
under the best seeing conditions (FWHM 0"55 . . .0."8). See 
also Fig.Q 
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Table 2. Positions and galactic reddening dSchlegel et al. 1 998) 
for the three COMBO-17 fields analysed. All observations 
were obtained at the Wide Field Imager at the MPG/ESO 
2.2 m-telescope at La Silla. 
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tral bins. This way, COMBO-17 provides identifications and 
reasonably accurate redshifts not only for galaxies but also 
for quasars, a novelty pioneered in CADIS (Wolf et al. 1999) 
and more recently applied in the SDSS ( Richa rds et al. 2 001 ) 
and to observations with superconducting tunnel junctions 
JdeBruiine et al. 20021) . 

All observations pres ented here were obtained with the 
Wide Field Imager ( WFI. iBaade et al. 1998b at the MPG/ESO 
2.2-m telescope on La Silla, Chile. The WFI provides a field 
of view of 34' x 33' on a CCD mosaic consisting of eight 2k 
x 4k CCDs with a scale of 0f'238/pixel. The total observing 
time for COMBO-17 was about ~175 ksec per field, including 
a ~20 ksec exposure in the R-band with seeing below 0f'8. 

Observations and data analysis have been completed on 
three fields including the Chandra Deep Field South (CDFS), 
covering an area of 0.78 deg 2 (Tab. 0. The deep /?-band im- 
ages have 5-0 point source limits of R ps 26 and provide the 
highest signal-to-noise ratio for object detection and position 
measurement among all data in the survey, except for L-stars 
and quasars at z > 5, w hich can be extremely fa int in R. 

Using SExtractor JBertin & Arnouts 1996I) . we obtained 
a catalogue of ~ 200,000 objects with positions and the 
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Fig. 1. COMBO-17 filter set: Total system efficiencies are shown in the COMBO-17 passbands, including two telescope mirrors, 
camera, CCD detector and average La Silla atmosphere. Combining all observations provides a low-resolution spectrum for all 
objects in the field. 



SExtractor photometry MAG -AUTO (for details on the con- 
struction of the catalogue see lWolf et al. 2001bl) . This .R-band 
selected source catalogue was then used to obtain SEDs across 
17 passbands with a photometric technique specifically tai- 
lored for measuring colour indices with high signal-to-noise 
and as little as possible interference by changing observing 
conditions. To this end, we projected the object coordinates 
into the frames of reference of each single exposure and mea- 
sured the object fluxes at the given locations using with the 
well-established MPIAPHOT approa ch of a se eing-adaptive 
weighted aperture jRoser & Meisenheimer 199lh . 

The COMBO-17 photometric calibration is based on a 
system of faint standard stars in the COMBO-17 fields, 
which we tied to spectrophotometric standard stars during 
photometric nights. Our stan dards were selected from the 
Hamburg/ESO survey database iWisotzki et al. 2 000 ) of digital 
objective prism spectra. By having standard stars within each 
survey exposure, we were independent from photometric con- 
ditions for imaging. More details of the data reduction will be 
provided in a forthcoming technical survey paper (Wolf et al., 
in prep). In the present paper, all magnitudes are quoted with 
reference to Vega as a zero point. 

3. The quasar catalogue 

3. 1. Sample definition 

The quasar sample is extracted from the full survey catalogue 
purely on the basis of spectral information, thus deliberately 
ignoring morphological evidence. We believe that the classifi- 
cation of the 'fuzzy spectra' based on 17 filters, outlined be- 
low, allows to clearly differentiate between stars, galaxies and 
quasars, providing a safer separation between the object classes 
than morphological criteria. This is particularly important in 
the context of our interest in low-luminosity AGN which may 
well appear extended on a deep 7?-band image. Our purely spec- 
troscopic classification approach ensures that the sample will 
not be heavily biased against such objects. 

The sample finally used for all analyses in this paper is de- 
fined by limits in magnitude and redshift. Objects are selected 
to have a magnitude of R > 17 to avoid the saturation regime 



in the individual frames, and R < 24 which is where the com- 
pleteness of the quasar identification has dropped below 30% 
(see Sect. I3.3I >. Since quasars show brightness variations, the 
sample selection will depend on the epoch of observation at 
the faint end. Here, we have used our R-band photometry from 
January 2000. 

The sample is further limited to redshifts of z > 1.2, since 
host galaxies may contribute significantly to the spectra at 
lower redshifts, where the 4000 A-break is still contained 
within the filterset. Our templates currently contain only pure 
quasar and pure galaxy spectra, but no mixed templates with 
contributions from both. Therefore, identification and redshift 
estimation of low-luminosity quasars at z < 1 .2 is not straight- 
forward at this point, and complicated completeness issues 
arise in the low-z domain. 

In order to ensure minimum contamination from non-active 
galaxies, we have set our probability threshold for an object to 
be classified as quasars quite high. As a result, we have elim- 
inated many trustworthy low-luminosity AGN from the sam- 
ple, but keep a well-controlled selection of higher-luminosity 
QSOs. We have chosen this conservative approach because we 
have not yet obtained spectroscopic confirmation redshifts for 
a sufficiently large sample to understand the selection function 
of low-luminosity Seyfert galaxies well enough. 

The catalogue resulting from this selection contains 192 
quasars between z = 1.2 and z = 4.8, with a median redshift 
of (z) = 2.23. Examples of quasar filter spectra are shown in 
Fig.|2j and the Hubble diagram for the full sample is presented 
in Fig. |3 

3.2. Classification and redshift estimation 

The photometric measurements from 17 filters provide low- 
resolution spectra for each object which are analysed by 
a statistical technique for classification and redshift estima- 
tion based on sp ectral template matching (for details see 
IWolfetal. 20"0l3) . Meanwhile, we have improved the tem- 
plate library for quasars by deriving it from the r ecent SDSS 
QSO template spectrum Jvanden Berk et al. 200 lh . rather than 
from the earlier used pure emission line contour by Francis 
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Fig. 2. Filter spectra of example quasars: The three panels show 
quasars at different redshifts across the range addressed in this 
paper. The filter spectra are plotted with horizontal bars resem- 
bling the filter width and vertical bars for la flux errors. Due 
to variability repeated observations in B and R filters can show 
different flux levels for the same object. Variability is compen- 
sated for the construction of SEDs using R-band images which 
are available for every observing run. The best-fitting template 
spectra at indicated redshifts are plotted as grey lines. 



et al. (1991). This leads to a better detection and redshift es- 
timation of quasars at z < 2.5 where the 'little blue bump' 
may render the spectral shape between the prominent emis- 
sion lines as quite different from the power-law which we as- 
sumed for the quasar templates previously. Since the CADIS 
work rfWolf et al. 19 99). we have also learned that photometric 
redshifts for quasars are strongly influenced by flux variability. 

Since the multi-colour observations were collected over a 
period of two years, and quasars as well as some stars show 
variability, it was necessary to correct for variability when con- 
structing the 17-filter spectra for classification. For this purpose 
we used the 7?-band observations which are available for each 
observing epoch, allowing us to completely map the variabil- 
ity at least in this band. When constructing the SEDs of vari- 
able objects, we related the measurements of observed pho- 
tometric band to the 7?-band magnitude obtained in the same 
observing run. As a result, the SEDs are not distorted by long- 
term magnitude changes. This variability correction works well 
only as long as flux variations are not accompanied by major 




Fig. 3. The present sample of 192 quasars: Distribution of red- 
shifts over observed R-band magnitude. Notice that the sample 
has been truncated at z — 1 .2 and R = 24. 



changes in spectral shape. Furthermore, we also cannot account 
for short-term variability on a time scale of days, as within each 
observing run we have obtained only one /?-band image. 

3.3. Completeness correction 

The subsequent analysis includes only objects with success- 
ful z estimates. It is therefore critical to understand for which 
quasars the data permit such a classification. Given the pho- 
tometric properties of the survey, we can produce mock cata- 
logues of stars, galaxies and quasars with realistic photomet- 
ric errors, and investigate the classification performance as a 
function of object type, redshift and magnitude. The 17 -filter- 
spectra in the mock catalogues are constructed from the li- 
brary templates, using the empirically derived observational 
error distributions. With this approach, the completeness of 
the classification and redshift estimat ion can be derived from 
Monte-Carlo simulations (see IWolf et al. 2001al) . We implic- 
itly assume that the spectral templates truly resemble observed 
objects, and if that is not the case, our simulation will be too 
optimistic. 

In fact, we can test whether our completeness maps are real- 
istic, given that within our fields and selection limits 12 broad- 
line AGN are known from spectroscopic observations in the 
CDFS (Hasinger 2002, priv. comm.; Szokoly et al, in prep.). 
The map predicts that 10 out of 12 objects should be identified, 
while in fact we recover 8. Two objects are missing from our 
sample although they lie in regions of high expected complete- 
ness: A QSO with broad absorption lines (BALs) at z = 3.6 and 
a Sy-1 galaxy at z = 1 .2. The BAL QSO was misclassified as a 
star because of its unusually star-like colours. The Sy-1 galaxy 
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resides just at our low redshift limit which we adopted to avoid 
incompleteness arising from the pure AGN spectra getting con- 
taminated by host galaxy contributions. We conclude that on 
the whole our map is realistic, but that (a) rare QSOs with un- 
usual colours could still escape our attention, and (b) right at 
the low-redshift limit the level of completeness might be re- 
duced compared to our maps. 

The product of the simulations for quasars is a complete- 
ness map (selection function), providing a formal probability 
C(R,z,SED) that a quasar of given intrinsic properties is re- 
covered, in bins of observed /?-band magnitude and redshift. 
Our quasar template library contains a range of emission line 
strengths and a range in spectral indices or continuum colours. 
Effective spectral indices depend on redshift as the quasar con- 
tinuum is not a pure power law. However, for a quasar at 
z = 2.0 the range of template B — I colours runs from about 
+0.35 to +1.75, corresponding roughly to power law indices 
of a = — 1 .66 . . . + 0.4. After we recognized that the simulated 
completeness depends very little on the spectral index, we col- 
lapsed the completeness function into C(R,z) removing the ex- 
plicit SED dependence (see Fig.|4j. 

At most redshifts, our classification algorithm is more than 
90% complete for quasars with R < 23. The contour line for 
50% completeness ranges around R ~ 23.5. The redshift de- 
pendence of the completeness shows conspicuous oscillations 
above redshift 2, which are caused by the strong signature of 
the Lyman-oc emission line migrating through the filter set and 
alternating between visibility in a medium-band filter and in- 
visibility when it falls between two neighboring medium-band 
filters. Whenever the line is visible quasars can be identified to 
fainter levels due to its contrast. 

At certain redshifts, the selection function reaches values 
above 100%. This occurs when redshift aliasing creates local 
overdensities in the estimated z-distribution based on a flat un- 
derlying simulated distribution. Such 'overcomplete' regions 
are always accompanied by adjacent 'undercomplete' zones, 
so that the total number of objects is conserved. 

3.4. Influence of limited redshift accuracy 

For the analysis it is crucial to check to which extent the limited 
redshift accuracy of a z « 0.03 (as determined from the sim- 
ulations) could affect inferences about luminosity functions. 
Two major aspects need to be explored, redshift aliasing and 
catastrophic mistakes, which are both irrelevant for interpret- 
ing well-exposed data from an aperture spectrograph, but could 
both play a role in our case: 

1. Aliasing results from structures which are finer than the 
redshift resolution, violating the sampling theorem; it ap- 
pears as fake structure on a scale that is typically slightly 
larger than the resolution. This effect is doubtlessly present 
in our data, as is also visible from the z-dependent oscilla- 
tions in Fig.|4] In this paper, we avoid dealing with the prob- 
lem altogether by only considering redshift bins of typically 
Az = 0.6, an order of magnitude larger than the resolution 
limit. 



2. Catastrophic mistakes occur in certain regions of colour 
space where different interpretations can be assigned to the 
same colour vector and probabilistic assumptions are used 
to make a final redshift assignment. Obviously, a number 
of cases could lead to the assignment of a wrong redshift, 
but their effect on luminosity functions should be relatively 
small unless a major fraction of objects were affected. The 
only real impact would be in a situation in which many 
low-redshift objects of medium luminosity were wrongly 
assumed to reside at high redshift, boosting the abundance 
of the rare luminous quasars. We see no evidence for any 
such effect to be important. 

A large-scale performance check of our spectrophotomet- 
ric redshifts based on a substantial sample of spectroscopic 
redshifts is still pending. However, we have been able to test 
our classifications and redshifts against two samples of albeit 
limited sizes that have spectroscopic data. In the Sll field, 
the COMBO-17 sample overlap s at its bright end with the 
'2QZ' quasar survey llCroom et al. 20011) . and in the CDFS 
we found that a number of spectroscopically identified X-ray 
sources coincide with faint COMBO-17 AGN (Hasinger 2002, 
priv.comm.; Szokoly et al, in prep.). Fig. [5] shows a compari- 
son of spectroscopic and multi-colour redshifts for all 22 ob- 
jects from the COMBO-17 QSO sample with R < 24 for which 
spectra exist. Only one of these 22 objects, classified as a QSO 
at z ~ 0.6 by COMBO-17 (i.e. outside the redshift range con- 
sidered in this paper), was actually a misclassification; a 2QZ 
spectrum identified this object correctly as a White Dwarf/M 
dwarf -pair. The other 21 objects were found to be QSOs, and 
their redshifts agree all within ±0.05. This comparison sug- 
gests that the accuracy of COMBO-17 redshifts is certainly no 
worse than the assumed a z = 0.03, and the rate of catastrophic 
mistakes is probably well below 10 %. 

The CDFS contains several further faint Seyfert galaxies, 
particularly at z < 1 , most of which were recently identified as 
optical counterparts to faint X-ray sources. These are very chal- 
lenging objects for the multi-colour redshift estimation, as their 
host galaxies will contribute significantly to their overall SED. 
Clearly, more work is needed to tackle these objects, which are 
therefore explicitly excluded from the scope of this paper. 

4. Number counts 

From the AGN sample and the completeness map described 
above it is straightforward to compute surface densities as a 
function of apparent magnitude, using the relation 

where the summation runs over all AGN brighter than R. Aq 
is the formal total survey area of 0.78 deg 2 ; multiplying this 
with the value of the completeness map C(Rj,Zi) for object i 
yields the estimated 'effective survey area' A,- for that particu- 
lar source. In other words, the decreased probability of detect- 
ing very faint AGN, particularly relevant for certain redshift 
ranges, is simply interpreted as being equivalent to a reduced 
survey area for such objects. 
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Fig. 5. Spectroscopic vs. multi-colour redshifts: 22 QSOs from the COMBO-17 sample at R < 24 have spectroscopic identifica- 
tions. 21 of them were found to be QSOs at redshifts within ±0.05 of the multi-colour estimate, one of them is a White Dwarf/M 
dwarf -pair estimated at z < 1 . 
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Fig. 6. Cumulative surface density of AGN as a function of R 
band magnitude, corrected for Galactic extinction. 

Results are shown in Fig. [6] and listed in Table [5] Besides 
the full sample of all z > 1 .2 AGN, we also show the number 
counts for subsamples split at z < 2.2 and z > 2.2, the typical 
high-redshift limit of UV excess surveys. The two subsamples 
show very different trends towards the faint limit: While the 
low-z objects dominate at brighter magnitudes, their surface 
density increases only slowly with decreasing flux level. On 
the other hand, the cumulative number of AGN with z > 2.2 is 
a strong function of magnitude, and quite well described by a 
power law 9l(<R) oc /? -0 - 31 , corresponding to a power law in- 
dex of — 1 .75 for the differential number-flux relation. Towards 



Table 3. Tabulated AGN number counts as a function of R band 
magnitude. The 2nd and 4th column give the actual numbers of 
objects in the current COMBO-17 sample (within 0.78 deg 2 ), 
while the 3rd and 5th column give incompleteness-corrected 
surface densities per deg 2 . 
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the faint end it even seems as if the slope might be increasing, 
but this is not formally significant. Nevertheless, for R > 23 it is 
clear that high-redshift AGN with z > 2.2 become the dominant 
population. 

5. Estimation of rest frame luminosities 

For each object in the sample we have individual SEDs from 
the 17-filter spectrophotometry at our disposal. In order to de- 
rive absolute magnitudes for the subsequent analysis, we de- 
cided to tie all luminosities to the UV continuum level at 
Arest = 145 nm. This was realised by integrating the best-fitting 
redshifted template over a synthetic narrow rectangular pass- 
band at 143 nm-147 nm, thus avoiding the Si IV/O IV emission 
line at 140 nm. This procedure enabled us to directly measure 
rest-frame luminosities over a redshift range 1 .4 <, z ^ 5 with- 
out any need for extrapolation. 



8 



Wolf et al.: Evolution of faint AGN from COMBO-17 



— i — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — r— 

1.5 '- - 

BO - 
03 

1 ( : . 

I 0.5 -_ — ; /■.'V\s:^/ ■ ; . ■ ■ ■ — -_ 

o Z . ' Z 

- 

« : . ■ " : 

— i i i i i i i i i i i i i i i i i i i i i 

1 2 3 4 5 

Redshift z 

Fig. 7. Derivation of the internal K correction used in this pa- 
per: Each point corresponds to one individual AGN measure- 
ment based on the SED fits within COMBO-17. The horizontal 
line is a linear fit, allowing to predict luminosities at 145 nm 
with an accuracy of 0.24 mag. 



N 



5 


1 i i i | i i / i | i i i 

/ 

/ 

/ 

/ 

/ ■ 


1 1 1 | 1 


4 

3 


/ 

/ 
/. 

/ ' - ". 
/' ' 

/-.•■'.' 




2 


h. ■-. ■„ " 

/ . - . 

/■;.;' ■ ... . . 

/ . ■ 
- / . 1 ■ • ■ 

/ . ' ■ • ! 

7 "■"'.*.■ 




1 


1 1 1 1 1 1 1 1 1 1 1 1 


~ 



-22 -24 -26 -28 -30 
M 145 [mag] 

Fig. 8. Distribution of the input sample over absolute magni- 
tudes and redshifts: The dashed lines indicate the imposed sam- 
ple limits of z > 1 .2 and R < 24. 



While such individually determined luminosities may suf- 
fer less from biases arising in uncertain assumptions on AGN 
spectral energy distributions than some previous samples, they 
considerably complicate the statistical exploitation. Although 
most of our objects are detected in all or nearly all of the 17 
filter bands, the sample is presently defined only in the 7?-band: 
(1) The flux limit is homogeneously truncated at R < 24. (2) 
The completeness map has the R band magnitude as one of 
its two independent parameters. It is also assumed in the algo- 
rithms used for luminosity function estimation that the survey 
is defined by just one flux limit. We have therefore explored 
whether it is possible to estimate absolute UV magnitudes M145 



from the measured R band flux, using a modified but basically 
traditional K correction approach. 

For this purpose we computed for each quasar its abso- 
lute magnitude directly from the distance modulus, Mo = R — 
(m — M), equivalent to assuming a K correction for a power- 
law spectrum with slope a = — 1 . We then plotted the differ- 
ences AM between Mo and M145 (the spectrophotometric esti- 
mates) against source redshift; this is shown as the distribution 
of points in Fig. This plot shows that at least at redshifts 
z > 2 there is a well-defined relation between z and AM. When 
fitting a constant line through these points, as shown in Fig.0 
we get an overall rms scatter of 0.24 mag (0. 1 3 mag for z > 2.2 
and 0.26 mag for z < 2.2), which can only insignificantly be 
improved by fitting e.g. a low-order polynomial. This line de- 
fines then an empirical correction relation K(z), so that for the 
redshift range of interest in this paper we can use R band mag- 
nitudes to estimate M145 with an accuracy of 0.24 mag, even 
without using any SED information. 

Most earlier studies of the optical AGN luminosity func- 
tion, especially those focussing on lower redshifts, have ex- 
pressed their results in terms of blue magnitudes Mb- In order 
to facilitate a statistical comparison, we obtained a crude esti- 
mate for the offset Mb — Mi 45 by repeating the above procedure 
for the rest-frame B band. Note that for all z > 1 .2 objects this 
involves a good deal of extrapolation, and the diagram corre- 
sponding to Fig.0shows more than 0.6 mag of scatter around 
the mean trend. Nevertheless, the offset of 1.75 mag thus de- 
termined is actually identical (to 0.01 mag) to the result when 
the same quantity is measured in the mean quasar energy dis- 
tribution of Elvis et al. (1994). In conclusion this means that 
all quoted absolute magnitudes M145 can be converted into Mb 
using the relation Mb = M145 + 1 .75. 

In the subsequent analysis we use two sets of cosmo- 
logical parameters, an Einstein-de Sitter universe with Ho = 
50 km s -1 Mpc~\ £l m = 1.0 and £l A = 0, and a flat 'con- 
cordance model' with Ho = 65 km s Mpc~', £2„, = 0.3 and 
£2a = 0.7. While the former is now physically almost obsolete, 
it still has some relevance as a 'yardstick' model, as most ear- 
lier studies of AGN evolution were expressed preferentially in 
these terms. Actually we found that the results change very lit- 
tle when switching between the two models, except for small 
shifts in both axes, and our displayed results always refer to 
the 'concordance universe' unless explicitly stated otherwise. 
Figure|8]shows the distribution of absolute magnitudes vs. red- 
shift for the entire AGN sample. In terms of the conventional 
(if arbitrary) distinction between high-luminosity quasars and 
low-luminosity Seyferts around Mb — —23, our sample covers 
just the region close to this dividing line. It is therefore the first 
optically selected AGN sample of substantial size to probe this 
important regime of nuclear activity in galaxies. It also matches 
very well the luminosity range cov ered by low-z AGN surveys 
that provide a z — reference (e.g.. Ko hler et al. 19971) . 
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Fig. 9. (a) Binned differential luminosity functions for six non-overlapping redshift shells. Only luminosity bins completely 
covered by the sample are shown. In addition to the data with Poissonian error bars, each panel features the 1.8 < z < 2.4 
luminosity function as reference, (b) Cumulative luminosity functions for the same redshift shells. 
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Table 4. Tabulated binned AGN luminosity function. 
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6. Evolution of the AGN luminosity function 

6. 1. Non-parametric estimates 
6.1 .1 . Method of computation 

We employ the usual 1 /V ma x estimator JSchmidt 19681) to 
give the space density contributions of individual objects. 
Luminosity functions are then readily obtained by forming the 
appropriate sums: Denoting the luminosity-binned differential 
LF as (j)(M) and the cumulative LF as \|/(M) we have 

M+AM z+Az i 

M— AM z-Az V > 
M z+Az i 

-°°z-Az V ' 

where each value of V, is the total comoving volume within 
which object ; would still be included in the sample. The value 
of Vj follows from integrating the survey selection function 
C(R,z), introduced in Sect. 13.31 over all relevant redshifts: 

/•Zmax dV 

V i (R,z)=A J C{R,z)—dz (3) 

fcf. lWisotzki 1 998) where Aq is again the formal total survey 
area, Aq = 0.78 deg 2 in the present case. The assigned error 
bars are based purely on Poissonian shot noise due to the lim- 
ited AGN counts within each bin, i.e. we ignored errors in mag- 
nitude or redshift estimation. The error bars for a given sum are 
then given by 

<^ = j£l/V?(M,z). (4) 

It is well known that binned estimates of luminosity functions 
suffer from a number of biases. In particular, two effects are 
worth being recalled, both of which can be seen as two dif- 
ferent variants of Malmquist bias and lead to the observed LF 
appearing flatter than the intrinsic one: 

- Binning in magnitude and a steep LF create a positive bias 
in the derived space densities, due to a shift in effective bin 
centres towards lower luminosities. 

- Differential evolution within redshift shells leads to an 
overestimation of space densities where the LF is steep. 

While the first effect is absent in the cumulative LF representa- 
tion, the latter is unavoidable as long as no model assumptions 



for the evolution of space densities are made. It should thus be 
kept in mind that the nonparametric LFs discussed in the fol- 
lowing section do not necessarily show the data 'as they are', 
but just form the most conventional way of presenting luminos- 
ity function data. 

6.1.2. Results 

Inside the redshift interval 1.2 < z < 4.8 covered by the 
COMBO-17 sample we defined six redshift shells of equal 
sizes to trace the evolution of the QLF Figure [9] gives a syn- 
opsis of the results; for the benefit of interested readers who 
wish to use our data for their own computations we present the 
tabulated binned LFs in Table [4] This table shows also how 
many objects contributed to each bin. The numbers given in 
this table show that the differences between adjacent redshift 
shells are not huge, i.e. evolution is relatively moderate over 
this redshift range. We therefore have deviated from the usual 
practice to plot all LFs into one frame, as any real trend would 
be very hard to discern in such a diagram. Instead, we created 
individual subpanels for each redshift shell and plot each lumi- 
nosity function separately, but together with the 1.8 < z < 2.4 
LF taken as a reference. This visualisation shows that the mea- 
sured values of the LF are below the corresponding reference 
values in nearly all data points. In other words, we detect an 
unambiguous maximum in the comoving AGN space density 
near these redshifts, and a significant drop towards both lower 
and higher z- The existence of such a maximum becomes even 
clearer in Fig. where integrated space density at given lower 
bound in luminosity (taken directly from the cumulative LFs in 
each redshift shell) is plotted directly against z- We return to 
a more detailed discussion of the evolution of space densities 
and LF shape properties below. 

6.2. Parametric analysis of the luminosity function 
6.2.1. Ansatz 

Within certain bounds and accuracy limits, an observed lumi- 
nosity function can usually be described quite well by some 
simple analytic expression, allowing one to compress the re- 
sults into a few well-determined numbers, with the positive 
side effect that the above mentioned Malmquist-type biases due 
to data binning can be avoided by obtaining luminosity func- 
tion and evolution parameters simultaneously from an observed 
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sample. This was first demonstrated by Marshall et al. (1983) 
for a very simple power law model of the evolving QLF. Since 
the parametric forms employed for the analysis in this paper 
are somewhat non-standard, we use the following paragraphs 
to spell out the adopted ansatz. 

Shape of the luminosity function: The most common ana- 
lytic description of the QLF is a double power law, involving a 
bright-end slope yi, a faint-end slope 72, and a smooth turnover 
at a characteristic 'break' luminosity M*. We have explored 
this ansatz and decided not to use it for the present sample, 
mainly because we fail to identify a well-defined break lumi- 
nosity in the data. Direct fitting of a double power law LF to 
any of the redshift shells subsamples shows that M* is a very 
ill-constrained quantity that often takes a value outside the lu- 
minosity range covered by the sample; furthermore, even mi- 
nor modifications in the subsample definition can cause major 
changes in M*. Since the luminosity functions in Fig. [9] un- 
doubtedly show some indication of curvature, we decided to 
parametrise the LF, at given z, as a polynomial in M145: 

log<|> = £A^. (5) 

1=0 

where /j = M145 — M* and M* is a constant which can be cho- 
sen freely but may also contain the effects of luminosity evo- 
lution (see below). In this representation, Aq = log(j)o is the 
space density normalisation at M145 = M*, which in our case 
has mostly been set to M* = —25 (roughly corresponding to 
Mb = —23.25). For n = 1, the LF is a single power law with 
slope y= -2.5Ai - 1. 

Evolution parameters: The last decade has seen a some- 
times hot debate over the question whether 'pure luminosity' 
or 'pure density' evolution were more appropriate. We have 
implemented both parametrisations as well as mixed forms. As 
independent variable we use 

C - log (6) 

which vanishes for z — Z re f • F° r our analysis we avoid extrap- 
olating to redshifts outside the sample definition interval and 
adopt z K f = 2. Pure density evolution is then expressed as poly- 
nomial expansion of the space density normalisation: 

log(|>o(z) = £c,C y . (7) 

;=o 

For mo = 1, the parameter C\ is equal to the density evolution 
index ko used in several earlier analyses. 

On the other hand, pure luminosity evolution involves a re- 
definition of the LF parameter M*: 

mi 

M*(z)=M k & = 0)+£B J y. (8) 

7=1 

(Note that unlike for a 'broken power law', M*(£ = 0) is here 
an arbitrary constant and not a free fit parameter.) Again, the 



case of mi = 1 corresponds to a standard form already used 
by several previous authors, M*(z) « (1 +z) kt - In this case the 
luminosity evolution index is Icl = 0.4Z?i. 

Fitting method : We follo wed a maximum likelihood ap- 
proach ( Mars hall et al. 19831) to search for optimal parameter 
combinations, using a modifie d downhill-simplex minimisa- 
tion scheme (Pre ss et al. 19921) . Several goodness-of-fit tests 
were applied to check whether the best-fit model was ade- 
quate. These included one- and two-dimensional Kolmogorov- 
Smirnov tests for the distributions over the (M,z) plane, as 
well as a % 2 test that directly compares observed and predicted 
numbers of AGN in redshift/luminosity bins. A fit was con- 
sidered satisfactory only when it passed all of the statistical 
tests with a probability for the validity of the null hypothesis 
of > 10%. This highly conservative approach ensures that the 
missing distribution-free properties of the KS test (see Wisotzki 
1998 for discussion on lLilliefors et al. 796 8) can be neglected. 

6.2.2. Results 

Starting with the simplest possible luminosity function, we 
found that a single power law LF as an overall shape model 
had to be rejected as there is weak but significant curvature 
in the LF, especially for the lower redshift ranges, z < 2.4. But 
already a 2nd order poynomial gives a statistically adequate de- 
scription of the LF at all redshifts. Note that this form involves 
one free parameter less than the case of a double power law LF. 

As a first evolution mode we explored pure density evo- 
lution (PDE). A global PDE fit over the entire redshift range 
1.2 < z < 4.8 was achieved with again just a 2nd-order poly- 
nomial to reproduce the maximum of comoving space den- 
sity around z — 2. Alternatively, the sample can be described 
equally well by pure luminosity evolution (PLE), albeit requir- 
ing one more parameter than PDE. Our best-fit PLE model fea- 
tures a 3rd order polynomial for the function M*(z). The pre- 
dicted parametric luminosity functions are compared with the 
binned nonparametric estimates in Fig. EH Both fits provide 
excellent descriptions of the data, with goodness-of-fit proba- 
bilities of p > 50 % for all statistical tests (see Tab.[5]for best- 
fit models). The differences between PDE and PLE luminosity 
functions become significant only outside the luminosity range 
sampled by the COMBO-17 data. 

In Fig. ^3 we show the cumulative luminosity function as 
a function of redshift, for different limiting luminosities, again 
together with the corresponding nonparametric estimates. For 
simplicity, only the PDE result is plotted. Since the fit is 
presently unconstrained below z — 1.2, we have marked the 
extrapolation into this region by a dotted line. The fit is com- 
pletely consistent with the individual binned datapoints. 

7. Discussion 

7. 1. The maximum of AGN activity 

The dominant feature in Fig.^^is the peak of comoving AGN 
space densities around z — 2. Although its existence was be- 
yond doubt for a long time, rarely has it been possible to locate 
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Fig. 10. Best-fit parametric representations of the evolving luminosity function, plotted against the binned LF data from Fig.|9ja). 
The density evolution model is shown by the solid, the luminosity evolution model by the dashed lines. 



Table 5. Coefficients of the best-fit analytic models describing the luminosity function and its evolution as either PDE or PLE. 
For details on the notation see text. 



Model Aq Ai A 2 M* 45 (0) C\ C 2 B { B 2 B 3 

PDE -5.600 0.2221 -0.03536 -25 0.3599 -15.574 

PLE -5.620 0.1845 -0.02652 -25 1.3455 -80.845 127.32 



the maximum within a single survey. UV excess-selected sam- 
ples typically reach up to z = 2-2.5 and mainly trace the rise 
at low z; conversely, most dedicated high-redshift surveys start 
being effective only beyond z — 3 and just show the decreas- 
ing branch. The COMBO-17 sample combines properties of 
both search techniques and covers enough redshift range that 
the peak of AGN activity is clearly bracketed. The PDE and 
PLE fits both place the maximum at a redshift of z max = 2.1. 

An important issue for the astrophysical interpretation of 
AGN evolution is a possible luminosity dependence of the peak 
location. For example, in hierarchical structure formation one 
expects objects of higher mass objects to form later, which in 
a simple scenario of correlated masses and luminosities should 
manifest in high-L AGN to show a peak at lower redshifts. This 
is certainly not the case in our data; on the contrary, one might 
be tempted to speculate from Fig.^Jthat there could be a grad- 
ual shift of the maximum towards higher redshifts when the 
luminosity limit is increased. 



We have tried to model such a trend by including explicit 
luminosity-dependent density evolution parameters, e.g. in the 
form of a linear dependency between absolute magnitude and 
Zmax- While the data are statistically consistent with a moder- 
ate shift, the best-fit model was always very close to simple 
PDE, despite the additional degrees of freedom. We conclude 
that there is no indication for a dramatic difference between the 
space density peaks of AGN of different luminosities, within 
the range covered by COMBO-17. 

A similar but much stronger effect, a shift of maximum 
space density towards lower redshift for very low luminosity 
AG N, has recently been reported for deep X- ray selected sam- 
ples (Co wie et al. 2 003: Hasinger et al. 20031) . While a detailed 
comparison between the properties of X-ray and optically se- 
lected AGN samples is beyond the scope of this paper, we just 
note two aspects which need to be taken into account: Firstly, 
the new X-ray selected samples probe even much deeper than 
COMBO-17 into the population of AGN with very low lumi- 
nosities. Assuming an Elvis et al. (1994) AGN spectrum, our 
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Fig. 11. Evolution of comoving AGN space density with redshift, for different lower luminosity limits, and for two cosmological 
models. Filled circles: Mus < —24; open circles: M145 < —25; filled squares: M145 < —26; open triangles: M145 < —27. The 
corresponding curves are integrated from the best-fit PDE models. 



low-luminosity limit of around M145 = —23 corresponds to a 
2-8 keV luminosity of roughly L x ~ 43.6; this is just the point 
where the difference becomes visible, according to Cowie et al. 
and Hasinger et al. 

But secondly, the X-ray and optical luminosity functions of 
high-redshift AGN also show significant dissimilarities. In par- 
ticular, the X-ray LF of type 1 (broad-line) AGN near logL t ~ 
44 is very nearly constant (cf. Cowie et al.; Hasinger et al.), 
while the optical LF around the corresponding M145 ~ —24.6 
is still significantly rising (cf. Fig. [5}. This may be indicative of 
different spectral shapes between low- and high-luminosity, or 
between low- and high-redshift AGN, all leading to nontrivial 
L x to M145 conversions. Thus, while on one hand our R < 24 
sample may still be not quite deep enough to see similar effects 
as found in the X-ray domain, it is not even obvious that exactly 
the same effects should be expected for a yet deeper sample. 

7.2. Evolution of the AGN luminosity density 

The contribution of AGN to the metagalactic radiation field is 
probably relevant for several wavebands. Recent satellite mis- 
sions were very successful in resolving the extragalactic X-ray 
back ground as ma inly due to distant low-luminosity AGN (e.g. 
iMivaii etaL 20001) . Likewise, the diffuse ionising UV back- 
ground is believed to be strongly influenced by AGN, but a 
quantitative synthesis still relies on several assumptions and 
extrapolations. One of the principal uncertainties has always 
been the shape of the low-luminosity end of the AGN luminos- 
ity function. In order to estimate the AGN luminosity density 
Pl(z), one has to evaluate the integral 

Pi(z) = jLi>(L,z)dL. (9) 




1 2 3 

redshift z 



Fig. 12. Integrated UV luminosity density Pi(i45) as a function 
of redshift, estimated from the cumulative LF (data points) and 
for both PDE (solid) and PLE (dashed) models. Below z = 1.2 
the relations are extrapolated (dotted lines). The discrepancy 
between PLE and PDE around the maximum is entirely due to 
the differences in the bright-end slope of the luminosity func- 
tion. 



The integral diverges only for a luminosity function slope 
y < —2, but the relative contributions of AGN at different lumi- 
nosity levels depend critically on the overall shape of the LF. 

The COMBO- 17 AGN sample goes deep enough that, for 
the first time, the quantity L(j) can be safely integrated without 
depending on heavy extrapolation into the unobserved range. 
We have computed both binned nonparametric as well as para- 
metric estimates of pi(z), which we show in Fig. ^1 We 
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present the monochromatic luminosity density Pi(i45), based 
on the value of XLx evaluated at A, = 145 nm (which in turn is 
directly derived fromM^s). Assuming a typical QSO spectrum 
such as the one given by Elvis et al. (1994), the quantity Pin^) 
can also be extrapolated into a frequency-integrated luminosity 
density. Since we are particularly interested in the AGN contri- 
bution to the cosmic production of hydrogen-ionising photons, 
we provide the necessary conversion: 

^<91.2nm = 1 -45 X A,Lx, = 145nm . (10) 

Notice that we have individual SED information on each of 
the COMBO-17 objects, which in principle could be used to 
derive a more accurate estimate of L; on . This will be done in 
a later paper dedicated to the synthesis of the UV background 
from the COMBO-17 AGN sample. 

For the nonparametric estimates we computed pi by sum- 
ming over the luminosity-weighted inverse volumes L/Vt up to 
the survey limit; the resulting numbers should set lower lim- 
its to the full integral over all luminosities. The two adopted 
parametric forms of PDE and PLE can be directly integrated 
to provide smooth functions Pl(z)- It is worth noting that the 
integral effectively converges within the COMBO-17 survey 
limits; the difference between setting the low-luminosity inte- 
gration boundary at M145 = —23 (approximately the COMBO- 
17 limit) and at M 145 = -10 (or L = 0) is only 2 %. In other 
words, the constraints on the faint end of the AGN luminosity 
function from COMBO-17 are already sufficient to estimate the 
total UV radiative output from AGN. This is also illustrated by 
the relatively good agreement in Fig. El between the binned 
estimates (which, as said above, are formally just lower limits 
to Pi) and the PDE model. 

On the other hand, the flat slope of the LF makes the con- 
tribution of brighter AGN to anything but negligible. In 
fact, the inclusion or exclusion of individual high-luminosity 
objects makes a noticeable difference for the binned estimate. 
(This is much less so for the parametric estimates, because 
of the uniform weighting inherent in the maximum likelihood 
procedure.) The importance of higher luminosity objects be- 
comes particularly apparent when comparing the PDE- and 
PLE-derived relations. The substantial and significant differ- 
ences seen in Fig.Elare not due to faint-end extrapolation ef- 
fects, but entirely originate in the much flatter LF slope of the 
PLE model at brighter magnitudes, already documented in Fig. 
E| above. We reiterate that these discrepancies are due to the 
limited power of the COMBO-17 dataset to discriminate be- 
tween evolution models. They will be resolved as soon as con- 
straints from other, brighter AGN surveys are combined with 
the COMBO-17 sample. 

7.3. Comparison with other surveys 
7.3.1. Redshift 1.5 <z<4.5 

In terms of the luminosity range, our quasar sample pushes 
to fainter limits than any previous survey. Our targets mainly 
have observed magnitudes of R = 20 . . . 24, while previous sur- 
veys either observed down to R SS 20 if they derived luminosity 
functions, or observed to R ^ 22.5 but produced object lists 



which did not constrain the luminosity functions very much 
further. Therefore, this work enters a new regime of studying 
low-luminosity AGN at intermediate to high redshifts. 

A first attempt to conduct an AGN survey with a strat- 
egy similar to COMBO-17 was performed, albeit on a much 
smaller scale, in the course of the CADIS survey (Wolf et al. 
1999), where a sample of 12 QSOs at z > 2 and R < 22 was 
identified within an area of 250 arcmin 2 . Their resulting sur- 
face density is twice as high as for COMBO-17, but because of 
the small sample size, the results of CADIS and COMBO-17 
are still formally compatible. The discrepancy underlines, how- 
ever, the importance of cosmic variance and hence the need to 
obtain large samples. 

Major previous wor k producing l uminosity functions in- 
clude the 2QZ at z ^ 2.3 JBovle et al. 20 00). the work by WHO 
at 2.0 < z < 4.5, and by SSG at z > 2.1 as well as the results 
from the high-redshift SDSS sample at z > 3.6 by Fan et al. 
(2001). As we have little overlap with all these previous sur- 
veys in luminosity, we cannot test directly to what extent our 
and their luminosity functions coincide. However, at the bright 
end our LF is completely consistent with an extrapolation of 
the SDSS-based LF given by Fan et al., with similar slope and 
similar normalisation. SSG cover a slightly wider range of red- 
shifts down to z = 2.7, but there still is not much overlap, so as 
before we need to extrapolate. The SSG slope is steeper than 
that of Fan et al., so that the SSG prediction for our regime is 
even further above the COMBO-17 results. 

On the lower-redshift side, our results smoothly connect to 
the 2QZ results. At intermediate redshift again, WHO used a 
broken power-law to characterize the QLF, where they found 
a steep end slope of a = —1.67 and a faint end slope of 
(3 = —0.45, which is mostly constrained by the probably non- 
optimal assumption of a broken power-law. The bright end of 
our LF varies around 0.4 to 0.7 in their units and can be con- 
sidered consistent with the WHO faint-end. 

7.3.2. Quasars atz>4.5 

The redshift range above 4.5 is subject to several studies from 
dedicated surveys to find faint high-redshift quasars. At high- 
est luminosities the SDSS has delivered sa mples of six ob- 
jects in 180 deg 2 at z* < 20 JFan et al. 200 lh an d 29 objects in 
~ 700 deg 2 at i* < 20.5 jAndersen etal. 200 lb . Medium-deep 
surveys across cover smaller areas, among which the BTC40 
has so far reported two z > 4.5 QSOs at / < 2 1 .5 across 36 deg 2 
and still have a list of fainter candidates at / ^ 22 to follow 
up JMonier et al. 200 21 The Oxford-Dartmouth Thirty Degree 
Survey (ODT) also aims at finding z > 4.5 quasars down to 
/ < 22 on 30 deg 2 (Dalton, priv. comm.). Monier et al. (2002) 
demonstrate the consistency of their currently available small 
numbers with the SDSS result and predict on the basis of the 
SDSS-LF a surface density of 0.026 /deg 2 at z > 4.5 ... 5.0 and 
/ < 19.9, and much less at z > 5. 

COMBO-17 covers only 0.78 deg 2 currently, but reaches a 
bit deeper. An extrapolation of the SDSS luminosity function 
suggests we should find 2.0 quasars per deg 2 at / < 23 and z > 
4.5. Our observation of some curvature in the fainter domain of 
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Fig. 13. Best-fitting parametric models of quasar space density in comparison: Solid lines are constrained by respective surveys 
and dashed lines are extrapolations to fainter luminosities. Left: COMBO-17 with 2QZ and SDSS (Fan et al.), at luminosities of 
M145 < [—28 . . . — 24]. The apparent disagreement of COMBO-17 with2QZ at brightest luminosities should not be taken seriously 
as the curves show the best-fitting model and not actual LF data. Right: COMBO-17 with WHO and SSG, at luminosities of 

M U 5 < [-28,-27,-26]. 



the luminosity function reduces that prediction to s=s 1 .25 / deg", 
so one or two objects is the total number to be expected in the 
COMBO-17 dataset. 

The present sample contains three objects at z > 4.5, all 
in the Sll field (see Fig. ^] for filter spectra). They all have 
/ < 22, although our selection should be complete to / < 23. 
Basically, this observation is statistically consistent with an ex- 
trapolation of both the SDSS LF and our own. As we can not 
draw strong conclusions about the cosmic abundance of these 
objects from COMBO-17, they have been excluded from the 
LF results presented above. 

7.4. Predictions for future deep surveys 

Given a parametrised model for the evolving AGN luminos- 
ity function, it is straightforward to predict surface densities 
(numbers of AGN per unit solid angle) for other survey spec- 
ifications. In Table [6] we provide a set of such numbers, com- 
puted for grid of R band magnitudes and redshifts. We deliber- 
ately stretch these predictions to the limit of credibility, e.g. in 
the lowest (R < 25) row or in the rightmost (z > 5) column, in 
order to enable a direct comparison with possible future ded- 
icated ultra-deep or very high redshift surveys. It should be 
understood that some of these numbers involve a considerable 
amount of extrapolation outside the areas covered by our data. 
For the same reason we list at each grid point the predictions of 
both PDE and PLE models, hoping that the difference between 
these approximately bracket the level of uncertainty, especially 
in the extrapolation regions. 



At intermediate redshifts and flux levels, well sampled by 
COMBO-17 sources, the agreement of PDE and PLE is excel- 
lent, simply reflecting the fact that both are valid descriptions of 
the observed data. On the other hand, the two evolution modes 
predict substantially different numbers of faint high-redshift 
AGN, the PLE prediction exceeding the PDE prediction by a 
factor of several. This is only partly due to the uncertainties im- 
posed by the limited dataset. Figure ITOl shows why this effect 
is actually expected from the properties of the two chosen evo- 
lution modes: A luminosity function which fits the data well 
at z — 2, which then is displaced either vertically or horizon- 
tally to account for the substantial negative evolution beyond 
z > 3, will result in dramatically different space densities of 
low-luminosity AGN. 

8. Conclusions 

We have presented work on the evolution of the quasar lumi- 
nosity function which has novel aspects by addressing two of 
the main unresolved problems in quasar research: The loca- 
tion of the peak of quasar activity, and the evolution of low- 
luminosity objects, the bulk of the quasar population: 

- Using the medium-band approach of COMBO-17 it was 
possible to break the colour degeneracies between stars 
and quasars around z ~ 3 which had prevented earlier op- 
tical surveys from mapping out the shape of the turnover 
in quasar activity from a single, homogeneous survey. As a 
result, we have measured a broad maximum around z ~ 2. 
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Table 6. Predicted AGN surface densities, based on (partly extrapolated) COMBO-17 results. Each field gives the cumulative 
AGN number per deg 2 , brighter than R and with redshift greater than z, as derived from the two simple evolution models discussed 
in this paper. 
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Fig. 14. Filter spectra of the three z > 4.5 COMBO-17 quasars 
in the present sample. See Fig.|2]for interpretation of the sym- 
bols. 



- By reaching faint magnitudes of R < 24, we were prob- 
ing the quasar LF and measured its slope down to rather 
low luminosity. The constraints on the faint end of the LF 
are fully sufficient to estimate the total UV radiative out- 
put. Although the two investigated evolution modes predict 
drastically different numbers of very faint (R <; 25), very 
high redshift (z 5) AGN, these have not much effect on 
the luminosity density integral. The main contribution to 
the UV background comes indeed from intermediate lumi- 
nosity objects within M145 w [—25, —28]. 



However, a number of questions and details remain still 
open. Among these are: 



1. The analysis of COMBO-17 catalogues themselves has 
so far omitted redshifts z < 1.2, because low-luminosity 
Seyfert galaxies, also those of type 1, may remain unde- 
tected due to their complex SEDs consisting of a mix of 
contributions from the host galaxy and the active nucleus. 
It will be a technical issue for the COMBO-17 classification 
to improve on that in the near future. 

2. At z > 4.5 COMBO-17 reaches to sufficient depth, but does 
not cover sufficient area to make a significant contribution. 
Wider area surveys such as the BTC40, the ODT or a deep 
sub-area of SDSS are needed to identify a larger sample of 
faint high-redshift quasars. 

3. Our sample alone fits well to both of the simplest evolu- 
tionary models for the quasar LF, the PDE and PLE mod- 
els. Thus, further improvements in determining the over- 
all shape of the LF and understanding its evolution on a 
broader basis will require the combination of the COMBO- 
17 sample with other samples at higher luminosity and 
lower redshift (to be addressed by a forthcoming paper). 

4. While our survey is conducted in the rest-frame UV do- 
main, one could wish to estimate the bolometric luminosity 
density by applying a bolometric correction factor. Using 
the Elvis et al. (1994) template SED, the relation would be 
^boi = 4.5 x Z445. However, we have little reason to assume 
a priori that the quoted bolometric correction is universally 
applicable to all AGN including our faint sample, so at this 
stage we refrained from following this up further. For sim- 
ilar reasons we do not want to derive accretion rates here. 

5. Probably, our sample breaks new ground in terms of depth 
and redshift coverage, however, it is clear that the results 
apply only to the evolution of broad-line type-1 AGN. 

6. Deep X-ray surveys with the Chandra and XMM/Newton 
observatories have started to charter into new territory of 
AGN research. Comparing the first results of X-ray and 
deep optical samples such as COMBO-17, we see similar- 
ities but also a number of differences. Exploring the inter- 
play between X-ray and optical domain for AGN will be a 
highly rewarding exercise, which will lead to a more com- 
plete understanding of the cosmic evolution of accretion- 
powered galactic nuclei. 
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In summary, understanding the accretion history of black 
holes in the nuclei of galaxies will certainly require some more 
quite fundamental work. However, establishing the evolution 
of faint quasars in the turnover epoch is one step forward. 
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